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The r/ c 7 transition form factor is calculated within a perturbative ap- 
proach. For the 7/ c -meson, a wave function of the Bauer-Stech- Wirbel 
type is used where the two free parameters, namely the decay constant 
f„ c and the transverse size of the wave function, are related to the Fock 
state probability and the width for the two-photon decay T[rj c — > 77]. The 
Q 2 dependence of the 7^7 transition form factor is well determined. 
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1 Introduction 

In 1995 the CLEO collaboration has presented their preliminary data on pseu- 
doscalar meson-photon transition form factors (see Fig. [I]) at large momentum 
transfer Q 2 for the first time [p]]. Since then these form factors attracted the 
interest of many theoreticians and it can be said that the CLEO measure- 
ment has strongly stimulated the field of hard exclusive reactions. One of the 
exciting aspects of the 717 form factor is that it possesses a well-established 
asymptotic behavior namely F W1 — > \[2f 1 ,jQ 2 where / 7r (= 131 MeV) is 
the decay constant of the pion. At the upper end of the measured Q 2 range 
the CLEO data |1||J] only deviate by about 15% from that limiting value. 
Many theoretical papers are devoted to the explanation of that little differ- 
ence. The perhaps most important result of these analyses, as far as they 
are based upon perturbative approaches (see e.g. 0-0]), is the rather precise 
determination of the pion's light-cone wave function. It turns out that the 
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pion's distribution amplitude, i.e. its wave function integrated over transverse 
momentum, is close to the asymptotic formf^] (~ x(l — x)). This result has far- 
reaching consequences for the explanation of many hard exclusive reactions in 
which pions participate (see, for instance, [jSHlOfl). 




Fig. 1. Meson-Photon transition form factors in e + e collisions. 

The situation is more complicated for the other cases, the 777 and the 7/7 
form factors. One has to determine not only the corresponding wave functions 
but as well the decay constants and the SU(3)p octet-singlet mixing angle for 
pseudoscalars. With the help of a few plausible assumptions a determination 
of these quantities from the 777 and 7/7 transition form factors seems possible 
I- 

There is a fourth form factor of the same type, namely the 7? c 7 form factor 
which is neither experimentally nor theoretically known. Since a measurement 
of that form factor up to a momentum transfer of about 10 GeV 2 seems feasible 
11 1 , a theoretical analysis and prediction of it is desirable. The purpose of this 
paper is the presentation of such an analysis. In analogy to the ttj case |)],|6| we 
will employ a perturbative approach on the basis of a factorization of short- 
and long-distance physics 0. Observables are then described as convolutions 
of a so-called hard scattering amplitude to be calculated from perturbative 
QCD and universal (process-independent) hadronic light-cone wave functions, 
which embody soft non-perturbative physics. The wave functions are not cal- 
culable with sufficient degree of accuracy at present and one generally has to 
rely on more or less well motivated model assumptions. 

In the case of interest the mass of the charm quarks, the 1] C meson is com- 
posed of, already provides a large scale which allows the application of the 
perturbative approach even for zero virtuality of the probing photon, Q 2 — > 0, 
and, therefore, our analysis can be linked to the two-photon decay width 



2 Allowing for a second term in the expansion of the pion's distribution amplitude 
upon the the eigenfunctions of the evolution kernel in order to quantify possible 
deviations from the asymptotic form, we find, from the recent CLEO data Q and 
within the modified hard scattering approach, a value of 0.0 ± 0.1 for the expansion 
coefficient B2 at the scale fi = 1 GeV. 
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r[?7c —* 77]- The experimental information on the latter width provides a con- 
straint on the t] c wave function. The valence Fock state probability P cS of the 
r] c , which is expected to lie in the range 0.8 — 1.0, offers a second constraint 
on the wave function and, for the simple ansatz we will use, determines it 
completely. We will show that variation of P c5 over the expected range has 
only a very mild influence on the final result, and hence our prediction for 
the transition form factor as a function of Q 2 turns out to be practically 
model-independent in the region of experimental interest, where potential Q 2 
dependence from higher order QCD corrections can be neglected. 

The organization of this paper is as follows: First we discuss the perturba- 
tive approach to the 7/ c 7 transition form factor, including the leading order 
result for the hard scattering amplitude and our ansatz for the wave function 
(sect. U). In the following sect. |3| the two parameters that enter our wave func- 
tion are fixed by relating them to the Fock state probability and the width for 
the two-photon decay. We present our results and conclusions in sect. f|. 



2 The perturbative approach 

In analogy to the case of the ivy case [|||| we define the 7^7 transition form fac- 
tor as a convolution of a hard scattering amplitude T# and a non-perturbative 
(light-cone) wave function \1/ of the r7 c 's leading cc Fock state 

KM) = jdx Jj^ H?(x, k ± ) T H (x, kj., Q) . (1) 


Here denotes the transverse momentum of the c quark defined with respect 
to the meson's momentum and x is the usual momentum fraction carried by 
the c quark. In contrast to the 717 case we do not include a Sudakov factor in 
eq. (|I|) and therefore we are not forced to work in the transverse configuration 
space. The Sudakov factor which comprises higher order QCD corrections 
in next-to- leading- log approximation [Jl2"yi3| , can be ignored for two reasons: 
First, due to the large mass of the c quark the QCD corrections only produce 
soft divergences but no collinear ones, and hence, the characteristic double logs 
do not appear. Secondly, the Sudakov factor is only relevant in the endpoint 
regions (x — > or 1) where it provides strong suppressions of the perturbative 
contribution. Since, however, the rj c wave function is expected to be strongly 
peaked at x — x , with x = 1/2, and exponentially damped for x — > 0, 1 the 
endpoint regions are unimportant anyway. 

The hard scattering amplitude in leading order is calculated from the Feynman 
diagrams shown in Fig. |2|. With one photon being almost on-shell qf ~ and 



3 



X X 

Fig. 2. The leading order Feynman diagrams contributing to the hard scattering 
amplitude for the r] c j transition form factor. 

the virtuality of the second photon denoted as q\ = —Q 2 , this leads to (with 
x — (1 — x)) 



r.(», k ± , Q) = - Q2 - ^ + pT) K - k , + (, » x) + 0(a,) (2) 

where M^ c (= 2.98 GeV) is the mass of the rj c meson, p is the ratio of the 
charm quark mass (m c ) and the r] c mass, for which we will take the value 
p = 0.5. The charge of the charm quark in units of the elementary charge is 
denoted by e c . Due to the symmetry of the wave function = ^(x), the 
two graphs provide identical contributions. 

For the i] c wave function we use a form adapted from Bauer, Stech and 
Wirbel PJ, 



*(x,k ± ) = ^|0(a;)£(k ± ). (3) 

f Vc is the decay constant (corresponding to f w = 131 MeV) which plays the 
role of the configuration space wave function at the origin. 4>(x) is the quark 
distribution amplitude which is parameterized as 



IX 



a) xx exp 



-a 2 Ml 

He 



Xo) 



(4) 



The normalization constant N^a) is determined from the usual requirement 
J dx(j)(x) = 1. The distribution amplitude (f|) exhibits a pronounced maxi- 
mum at Xo and is exponentially damped in the endpoint regions. This feature 
of the distribution amplitude parallels the theoretically expected and exper- 
imentally confirmed behavior of heavy hadron fragmentation functions. Fur- 
thermore, S is a Gaussian shape function which takes into account the finite 
transverse size of the meson, 



S(k ± ) = 167r 2 a 2 exp[-a 2 k; 



d 2 k ± 

167T 3 



S(kj 



(5) 
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Frequently used and for light mesons even mandatory [|15| ,[7| is a form of the 
kj_ dependence like exp[— b 2 kj_/xx}. Due to the behavior of the distribution 
amplitude (||) any explicit appearance of x in E can be replaced by x to good 
approximation. 



3 Fixing the parameters 



Let us start with the determination of the 7] c decay constant, a parameter 
which is not accessible in a model-independent way at present. Usually, one 
estimates f Vc through a non-relativistic approach which provides a connection 
between f Vc and the well-determined decay constant of the J/ip. We note, that 
the non-relativistic approach, which is only valid for Q 2 <C M 2 , is consistent 
with our definition of f Vc (see (0-0)) if relativistic corrections are ignored. In the 
non-relativistic approach the partial widths T[rj c — ► 77] and T\J/ij) — > e + e~] 
are related to each other 



rw - 



77] 



3 eta 2 



e e 



m 2 



Rs(o)\'' 



a* 

-3.4 — 

7T 

a, 
5.3 — 

7T 



(6) 



Here Rs{r) is the common non-relativistic S*-wave function of the J/ip and r) c 
meson, Ai^ parameterize the leading relativistic corrections, and the a s cor- 
rections have been calculated in [IB|. The wave function at the origin Rs{0) 
is related to the decay constants, and in the limit v 2 — > 0, a s — > one has 

77]/r[J/V 



U = fj/i> = p/2m c Tx\R s m and T[rj c -> jj]/T[J/iP -> e+e"] - 3e 2 c . The 
latter decay constant is model- independently determined from the J/ip lep- 
tonic decay width 



47re c 2 a 2 / 2 



e e 



3M 



5.26 ± 0.37 keV [17] 



(7) 



which leads to fj/^ = 409 MeV. However, the a s corrections in (|5|) are large 
(depending on the value of a s one prefers), and the relativistic corrections are 
usually large and model-dependent (e.g. Chao et al. |TS[ find A x = 5/12, A 2 = 
11/12 from a Bethe-Salpeter model). Estimates of the corrections typically 
lead to |H fjfj/j, = 1.2 ± 0.1 and T[ Vc -> 77] = (5 - 7) keV. 



The parameters entering the wave function are further constrained by the Fock 
state probability 
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i>p cc -=/^i^,k ± )i^^.^u + 



16tt 3 1 v ' 3 v^F V a 2 M 2 



Vc 



As we said in the introduction, one expects 0.8 < P cS < 1 for a charmo- 
nium state (for smaller values of P c5 one would not understand the success 
of non-relativistic potential models for these states). Since the perturbative 
contribution to the r/ c 7 form factor only mildly depends on the value of P cd , 
as it will turn out below, we use P c5 = 0.8 as a constraint for the transverse 
size parameter a. For f Vc = 409 MeV this leads to a = 0.97 GeV -1 , a value 
that is consistent with estimates for the radius (r 2 ) = 3 a 2 ~ (0.4 fm) 2 or the 
quark velocity v 2 = 3/(Ma) 2 ~ 0.3 from potential models JHJ]. 

The two photon decay width T[i] c — > 77], the experimental value of which 



still suffers from large uncertainties ||17|| ) can be directly related to the 7/ c 7 
transition form factor at Q 2 = 



rr 1 ™ 2 K 1 f rmi 2 7 - 5± ^ keV (direct) m 

nVc -> 77] = — 1^7(0) I = i r tot (9) 

4 [(4.0±1.5keV)- T 3^ v (BR) 

One may use this decay rate as a normalization condition for F Vcl (Q 2 = 0) and 
present the result in the form F Vcl (Q 2 ) / F Vcl (0) . In this way the perturbative 
QCD corrections at Q 2 = to the ?7 C 7 transition form factor, which are known 
to be large (see eq. @), are automatically included, and also the uncertainties 
in the present knowledge of f Vc do not enter our predictions. 



4 Results and Conclusions 



Let us turn now to numerical estimates of the 7^7 transition form factor. 
The left hand side of Fig. ^| shows that form factor for two different values 
of the Fock state probability P c5 . As already mentioned, we observe that the 
dependence on P c5 is weak. On the right hand side of Fig. |3| we present the 
result for the transition form factor Q 2 F Vcl scaled to a partial width T[r] c — > 
77] of 6 keV. 

For the values of the meson mass and the transverse size parameter a that 
we are dealing with (i.e. (aM^J 2 = 9) it makes also sense to consider the 
peaking approximation in which the hard scattering amplitude is evaluated 
at the position of the maximum of the distribution amplitude. The peaking 
approximation is formally equivalent to the replacement of the distribution 
amplitude (||) by a 5 function at X — Xq (into which it collapses in the limit 
(aM Vc ) — > 00). This approximation is numerically quite reliable and allows 
one to discuss the qualitative features of the model results in a rather simple 
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r\ c y form factor r\ c y form factor 

(normalized to r {r\ c ->yi\ = 6 keV) 




5 10 5 10 15 20 

Q 2 [GeV 2 ] Q 2 [GeV 2 ] 

Fig. 3. Left: The results F Vcl (Q 2 ) / F Vcl (0) for different values of the Fock state 
probability P C c as well as the approximation (|i"o| ) (with (kj_) — > 0) with and without 
a s corrections. Right: The predictions for Q 2 F VcJ (Q 2 ) scaled to T[r] c — ► 77] = 6 keV 
in the leading order of the perturbative approach (for P q g=0.8). The dashes indicate 
the Q 2 region where QCD corrections may alter the predictions slightly. 

fashion. By means of the uncertainty principle (aM Vc ) 2 3> 1 can be turned into 
(k±) <C M v 2, and hence one may also neglect the k\ dependence in the hard 
scattering amplitude (collinear approximation). Including l/aM Vc corrections 
to both, the peaking approximation and the collinear approximation, and using 
for the mean transverse momentum the relation 1/a 2 = 2(kj_) (see @) one 
arrives at the approximate result for Q 2 < M 2 



)2,^ ^ e cfvc _ ^7(0) 



mW ] ~ Q 2 + M 2 c + 2 (ki) " 1 + Q 2 /(M 2 c + 2 (ki» • [W) 

which agrees with the perturbative result to order (1 / aM Vc ) 2 . Eq. (|T0|) reveals 
that, to a very good approximation, the predictions for the 7/ c 7 form factor are 
rather insensitive to the details of the wave function. Only the mean transverse 
momentum following from it is required. 

To assess the quality of the approximation (|i~0|) we compare it for the special 
case of (k 2 ^) = to the full result from the perturbative approach in Fig. ^ 
(left hand side). We observe that, with increasing Q 2 , the two results growingly 
deviate from each other, at Q 2 = 10 GeV 2 the difference amounts to 10%. If 
one uses our estimated value of (k±) in ( |T0D the deviation from the full result 
is further reduced and amounts only to 4% at Q 2 = 10 GeV 2 . This little 
difference is likely smaller than the expected experimental errors in a future 
measurement of the ^7 form factor (see [0). These considerations nicely 



illustrate that the Q 2 dependence of the 7^7 form factor is well determined. 
The main uncertainty of the prediction resides in the normalization, i.e. the 
Tj c decay constant or the value of the form factor at Q 2 = 0. 
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Eq. (|10D resembles the Brodsky-Lepage interpolation formula proposed for 
the 7T7 transition form factor as well as the prediction from the vector 



meson dominance model (VDM). Our value of \jM 2 c +2(k 2 ± ) is 3.15 GeV 
which is very close to the value of the J/ip mass that one would have in- 
serted in the VDM ansatz^]. In the VDM the rjcj form factor at Q 2 = 
is given by F^° M (0) = e c gj/^ fj/^/Mj^ where the J/^fjc^f coupling con- 
stant can be obtained from the radiative decay J/ip — > r] c j fl7|j2"3|l . One finds 
F™(0) = 0.048 MeV" 1 and hence T VDM [r] c -> 77] = 2.87 keV which appears 
to be somewhat small as compared to the experimental values quoted in (|9|). 
Inclusion of a similar contribution form the ip' pole does not improve the VDM 
result since the ip' contribution is very small. In the case of two virtual photons 
q\ ^ the perturbative prediction F(g 2 , gf ) oc l/(—qf — q\ + M 2 c + 2(kj_)) 
differs substantially from the VDM. 

Let us briefly discuss, how a s corrections may modify the leading order result 
for the rjc'-f form factor: One has to consider two distinct kinematic regions. 
First, if Q 2 < M 2 c one can neglect the evolution of the wave function, and one 
is left with the QCD corrections to the hard scattering amplitude Th, which 
have been calculated in the peaking and collinear approximation to order a s in 
[24|. For the scaled form factor the a s corrections at Q 2 and at Q 2 = cancel 
to a high degree, and even at Q 2 = 10 GeV 2 the effect of the a s corrections is 
less than 5% (see left side of Fig. |3[). 

Secondly, for Q 2 3> M 2 c , one can neglect the quark and meson masses and 
arrives at the same situation as for the pions. The a s corrections to the hard 
scattering amplitude and the evolution of the wave function with Q 2 are known 
p| j2^ , |25|| . For very large values of Q 2 the asymptotic behavior of the transition 



form factor is completely determined by QCD, since any meson distribution 
amplitude evolves into the asymptotic form <f>(x) — > (fias(x) = 6xx, 



f v AQ 2 ) - ^ / dx ^ -> |^ , (g 2 - 00) (ii) 


The value of the moment (x^ 1 ) = J dx<()(x)/x evolves from 2.5 to the asymp- 
totic value 3. We note that the asymptotic behavior of the peaking approxima- 
tion (|T0l) is 16/9 f Vc /Q 2 - The deviation from ([TT|) demonstrates the inaccuracy 
of the peaking approximation for broad distribution amplitudes. 

A precise measurement of the strength of the ^7 transition form factor may 
serve to determine the decay constant f Vc (see, e.g. ( [TOP ). Though attention 

3 Considering the uncertainty in the mean transverse momentum and in the c-quark 
mass (see (|2|)), we estimate the uncertainty in the effective pole position following 
from (|io|) to amount to about 5%. 
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must be paid to the fact that the obtained value of f Vc is subject to large 
QCD corrections (about of the order 10-15% for Q 2 < 10 GeV 2 ) which should 
be taken into account for an accurate extraction of the rj c decay constant. 
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